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Abstract 

We discuss a dynamical matrix model by which probability distribution is associated with Gaus- 
sian ensembles from random matrix theory. We interpret the matrix M as a Hamiltonian repre- 
senting interaction of a bosonic system with a single fermion. We show that a system of second- 
quantized fermions influences the ground state of the whole system by producing a gap between 
the highest occupied eigenvalue and the lowest unoccupied eigenvalue. 
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I. INTRODUCTION 



Random matrices [1] first studied by Wigner [2] have applications in many branches of 
(many-body) physics, e.g., nuclear and molecular physics. The matrix elements are consid- 
ered statistical variables and take random values. On the other hand, one can consider a 
dynamical matrix model in which the matrix elements are dynamical variables so that the 
whole matrix has become a mechanical system, which may even be quantized and corre- 
sponding quantum field theory may be formulated. Such dynamical matrices have found 
use in e.g. description of disordered system and in string theory [3]. 

In the present work we treat a dynamical matrix model by which probability distribution 
in the ground state is associated with Gaussian ensembles from random matrix theory. We 
interpret the matrix M as a Hamiltonian representing interaction of a bosonic system, such 
as for example the nuclei in a molecule, with a single fermion, e.g. a single electron in 
a molecule. In particular, an eigenvalue of the matrix correspond to an eigenenergy of a 
fermion. Without going into so much detail of the single fermion interaction, we consider 
a backreaction from the system of the second quantized fermions. A characteristic effect 
of such a backreaction is to drive the matrix in the direction of lowering the energy of the 
eigenstates occupied by fermions. In the case of a molecule this backreaction is the force 
from the electrons pushing on the nuclei seeking to drive them into such a position so as to 
lower the filled single electronic orbit energy eigenvalues. Since a similar push to lower the 
unoccupied levels is absent, or even has the opposite sign, as will be argued later, hereby 
can easily arose a gap between the highest occupied eigenvalue (homo, where mo stands for 
molecular orbit) and the lowest unoccupied eigenvalue (lumo). The gap is what is called 
the homolumo gap [4]. The interest in the homolumo gap steams from the fact that the 
details of the model seem very unimportant. You should get it whenever you have a system 
of "bosonic" variables interacting with the fermions, provided that the "bosonic" variables 
are sufficiently soft as to yield to the pressure from the fermions. 

The purpose of this paper is to evaluate how distribution of eigenvalues of the matrix 
adjust to produce a homolumo gap in quite general setting. First, we formulate field theory 
corresponding to the large N limit of the random (N x N) matrix. Then we introduce 
the dynamical matrix model by which probability distribution in the ground state coincides 
with the large N limit of the 'free' Gaussian ensemble from the random matrix theory. The 
notion of the dynamical model enables us to model the interaction potential which arose as 
a consequence of the backreaction of the second quantized fermions. Finally, exploring the 
effect of adding this interaction potential to the 'free' matrix model, we obtain characteristic 
homolumo gap. 

II. FIELD THEORY FORMULATION OF THE RANDOM MATRIX MODEL 

The field theory we discuss is defined by the following functional integral 

Z[p\ — J Dp exp|(A — 1) J dxp(x) In p(x) + A J dxdyp(x)\n(x — y)p(y) — J dxp(x)P(x)^ , 



2 



for some (polynomial) function P(x) and a dimensionless parameter A. The integral (1) 
corresponds to the large N limit of the matrix integral from random matrix theory 

Z = J dMe~ TrP( - M \ (2) 

where P(M) is a polynomial in N x N matrix M. In terms of eigenvalues Xi of the matrix 
M, the matrix integral (2) can be expressed as 

Z= fl[dx i ]l(x i -x j ) 2X e-^ Xi \ (3) 



i i<j 



Here, the parameter A = 1/2, 1, 2 for real- symmetric, hermitean and quaternionic-real matrix 
M, respectively. The origin of the first two terms in (1) is the new invariant measure resulting 
from changing the variables from Xi to p(x), see Refs.[5, 6]: 

II(«. - - exp { A / ^) in p W + A / dM H, - y)M } . (4) 



Kj 

N 



(5) 



While in (3), the statistical variables are the eigenvalues x i: in the large N limit one in- 
troduces the density field p(x) which then by itself becomes statistical variable. In parallel 
with Ref. [7], we interpret (1) as a quantum field theory describing statistical random matrix 
model (2). By varying (1) with respect to p and taking derivative with respect to x we find 
that the most probable configuration satisfies following equation 

(\-l)^-2\Trp H = d x P, (6) 

where f H denotes Hilbert transform of the function /. From (6) we obtain the Riccati 
differential equation for the new field W = p H + ip in analogy with [8]: 

(A - l)d x W - XirW 2 - Wd x P = (pd x P) H - p H d x P (7) 

This equation corresponds to the usual Riccati equation for the resolvent function in the 
random matrix theory except that here the first term vanishes exactly for the hermitean 
matrix model (A = 1). The difference comes from the additional term in measure (5). 

Before introducing the matrix model which describes interaction of fermions and bosons 
we need the notion of dynamical system to be able to define and model the potential, 
and also, to be able to interpret the integrand in (1) as a probability density functional in 
the ground state of this dynamical model. This requirement is satisfied by the following 
Hamiltonian (see [8] for details): 

H = ^J dxp(x)A\x)A(x) + E , (8) 

where A(x) is given by 

51n$ 

A(x) = d x n(x) + id xj^, (9) 
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for $ being the ground-state functional (i.e., the square root of the integrand in (1)) and 
7r(x) the canonical momentum satisfying [d x ir(x), p(y)] = —id x 5(x — y). Now, we interpret 
Eq.(6) as a BPS equation of motion. Let us consider a simple example of the 'free' Gaussian 
ensemble defined by (3) with P(x) = x 2 . The operator A(x) of the corresponding dynamical 
model is 

A{x) = d x 7r(x) + i (^^p " ^P H {x) - *) , (10) 
and the Hamiltonian can be written as 

H= l -J dxp(x)(d x n(x)) 2 + l J dxp(x) ( A ~ 1 dx ^ - Xnp H (x) - xj + E , (11) 

The second term is called effective potential V e s and represents the physical potential for 
the dynamical model in question. The ground-state functional satisfying A(x)$ = is 

$ = cxp|- — - — - J dxp(x) In p(x) + ^ J dxdyp(x)\n(x — y)p(y) — — J dxp(x)x 2 1 , (12) 

giving the integrand of the (1) as a probability density functional in the ground state. The 
semiclassical solution is given by the solution of the equation (6), which in A = 1 case reads 
—irp H (x) = x and gives the usual Wigner semicircle law for the distribution of eigenvalues 
for the hermitean Gaussian ensemble. In the rest of the paper we restrict ourselves to the 
A = 1 case. 



III. INTERACTION WITH FERMIONS 

It is simply the idea of the interaction of the already described dynamical matrix model 
with the system of fermions that we postulate that dynamical matrix itself is the Hamiltonian 
for a single fermion. Then on top of that we second- quantize the fermions so that it becomes 
possible to have filled or empty all the (basis) states for this matrix. If we especially seek the 
ground state of the whole system - dynamical matrix plus fermions - we shall be interested 
in the case when the eigenstates of the dynamical matrix are filled below a certain value xp, 
the fermisurface (value), while the ones above xp are empty. In case we could ignore that 
this fermisurface (energy) depends on the (dynamical) state of the matrix the interaction 
with the fermions simply gives an extra potential 

Vfermion(^) = xQ{x F - x) . (13) 

At this point we introduce the zero-point energy for fermion modes. If we want to write 
a hermitean Hamiltonian contribution from a single particle state as a bilinear expression 
in the annihilation a and creation operators and symmetrized in taking the product we 
must take the commutator expression 

H=^[a\a]=^a-^. (14) 

The extra term —u/2 compared to the pure number operator term, coN = ua^a, is analogous 
to the zero-point energy term ou/2 for a boson model 

H =^-{ a \a} = coa t a + ^-. (15) 
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With this analogy between bosons and fermions in mind we can claim that there is a zero- 
point energy — 00/ 2 for fermion mode meaning that the energy of an empty level is indeed 
—uj/2, rather than the naive zero. Then of course the energy of the filled level should be 
u/2. 

The system which we consider of a dynamical matrix M interacting with a system of 
second quantized fermions described by annihilation Oj and creation a - operators denoted 
by the same indices {i,j} as the columns and rows in the dynamical matrix M. This means 
- with the convention of inclusion of the zero-point energy - that the energy term for the 
interaction of the dynamical matrix M with the fermions, and that includes actually all 
energy of the fermions, becomes 



Him. fermion = ^ [(«!, 4, . . . , O^), M(t) 



/dl\ 

\a N J 



(16) 



If we let Oj and a\ be annihilation and creation operators for the eigenstate of matrix M 
instead of for the original basis vectors, Hamiltonian (16) would reduce to 

#int. fermion = ^ a], Ct;j = ^£;(aja; - ]-). (17) 



In the following it is understood that the energy eigenvalues have been enumerated in in- 
creasing order, so that x\ < x 2 < . . . < Xi-i < Xi < Xi+i < . . . < xn- We can of course 
consider any energy number xf obeying x n . < xf < x nj+ i the fermisurface (energy), rif be- 
ing the number fermions. Instead of keeping the number of fermions fixed to nj, we include 
a chemical potential term in the Hamiltonian, 

H ch = —xp ■ #fermions = — xp • nf — —x F a\di. (18) 

i 

With such a term one arranges the minimal energy situation to have precisely the (wanted) 
value xp for the Fermi surface. Analogously to the fermion interaction term (17), we shall 
also for this chemical potential term choose symmetrized expression [aj,aj] in the annihila- 
tion and creation operators. If we consider the situation in which the fermions for a given 
state of the dynamical matrix M had adjusted to lower the energy most possible (for fixed 
M) then we would have the rif lowest eigenvalues filled and remaining N — nj eigenenergy 
levels empty. Including zero-point energy the value of the H int f crm i on part of the Hamiltonian 
would be 

n/ 1 N 1 If 

H^t. fcrmion I minimal = -Xi - y~] -Xi -> -- / dxp(x) assign (x F - x) , (19) 

fcrmion ' * Z Z Z J 

energy j=l i=Uf + l 

and the chemical potential term would give 

N 



1 1 x F f 

-f^ch | minimal = ~ ~X F + ~X F — > — / d,Xp(x) Sign (x F - x) . 

fermion Z . Z 2 J 



(20) 



energy i=\ i=nj + l 
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IV. GROUND STATE AND THE SPECTRUM DISTRIBUTION 



The full Hamiltonian for the system would contain additional terms coming from fermion 
interaction and chemical potential terms 



-^full H ~t~ -^int. fermion | minimal ~\~ H c \± | minimal 

fermion fermion 

energy energy 



(21) 



Properly it should, however, be stressed that although we shall formally use this expression 
it is in fact only physically justified when one seeks the ground state. If the system gets 
excited we also expect the fermion part of it to get excited and then one should in principle 
treat the fermionic part of the system as a properly second-quantized system. 
In seeking the minimum of effective potential of our dynamical model 



eff 



7T 



</.(• { —p 3 (x) + uj 2 x 2 p(x) — fip(x) — (xp — x)sign(xp — x)p(x) 



(22) 



and solving 5V e s/Sp = one obtains the semiclassical ground-state eigenvalue density 

p(x) = —\/p+(xf — x)sign(x_F — x) — uo 2 x 2 . (23) 



For x < xf, i.e., for filled states we have 



UJ U + Xp 
p{x) = -t 



7T 



X + 



2lu 2 



(24) 



and for empty states, x > xp, we have 

UJ I fl — Xp 



p{x) 



— X 



71 
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2u 2 



(25) 



Pictorially, we obtained semicircle law for both, filled and empty states, separated by the 
gap, when the chemical potential // = // + (2u)~ 2 satisfies xp < jl < xp + 1/cu 2 , and 
< Xp < 1, see Fig.l. Notice that the term (x — XF)sign(x — xf) in V e s has a constant 



p(x) 




FIG. 1: On the left we have filled levels and on the right empty ones. 

derivative, corresponding to a constant 'electric field', on each side of the Fermi surface xp. 
This 'electric field' pulls the energy levels away from xp, producing the homolumo gap. If 
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we didn't symmetrized the interaction with respect to filled and empty states, we would 
have 6{x — xf) instead of sign(x — xf), i.e. we would still obtain the gap, only it wouldn't 
be symmetric with respect to xf- 

A more intuitive physical picture could be, however, obtained in the following way. As 
we are interested only in the ground state we approximate the fermionic interaction with a 
smooth polynomial expression, as depicted in Fig.2. Choosing f(x) = ax 6 — bx 4 — ex 2 we 

f (X) 




FIG. 2: Smoothing the potential: upper graph is f(x) = x 2 — \x\, and lower is f(x) = —x 2 —x 4 +2x 6 . 



obtain for the ground-state functional 



$ = exp {\ f dxdyp ^ ~~ y)p(y) ~\f dx p( x )(f x4 - gx 2 ) I , ( 26 ) 



This correspond to 



J i i<j 



Xj) 



(27) 



giving, for deep enough wells, the two-cut solution. This shows that a homolumo gap can 
appear in a very generic system. Maybe even too generic, as it is basically just two-cut 
solution of a hermitean matrix model. However, the physical interpretation and field-theory 
formulation that we presented enables one to go beyond the ground state using rather 
standard techniques. Remember that it has been shown [1] that there are a few eigenvalues 
left outside of the semicircle, meaning that there are few states left in the gap. Thinking of 
our model as a field-theoretical model one can look for the instanton contributions [9] which 
are expected to give strongly-reduced but non-zero level density near the Fermi surface. 
Detail analysis of these effects is left for further investigation. 

One purpose of the presented studies in a very general setting is the application of the 
homolumo gap effect in the project of Random Dynamics [10]. There, one starts from the 
observation that at the present the energies at our disposal are extremely low, compared 
to the supposed fundamental energy scale, presumably to be identified with the Planck 
scale. Thus the 'poor physicist' can only hope to bring say a fermion from just below the 
Fermi-energy to just above it from the fundamental scale point of view. Obviously, the 
appearance of a homolumo gap of the fundamental energy scale order of magnitude would 
prevent creation of any fermion whatsoever. What we are really interested in is obtaining 
strong reduction of level density near the Fermi surface, instead of a genuine homolumo 
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gap. Namely, if the level density in the energy range we consider is especially low, then 
whenever a potential scattering from one momentum eigenstate to one with a different eigen- 
momentum should take place there will be for it especially small phase space available. That 
is to say there would be very few states to scatter into and thus we expect much fewer such 
momentum violating scatterings would take place. This strong enhancement in the accuracy 
of momentum conservation would be interpreted as a Random Dynamics spirit derivation 
of momentum conservation. 
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